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*5 ' Abstract 

In addition to thermal noise, which is essential to promote conformational transitions in biopoly- 
mers, cellular environment is replete with a spectrum of athermal fluctuations that are produced 
from a plethora of active processes. To understand the effect of athermal noise on biological pro- 
cesses, we studied how a small oscillatory force affects the thermally induced folding and unfolding 
transition of an RNA hairpin, whose response to constant tension had been investigated exten- 
sively in both theory and experiments. Strikingly, our molecular simulations performed under 
overdamped condition show that even at a high (low) tension that renders the hairpin (un)folding 



improbable, a weak external oscillatory force at a certain frequency can synchronously enhance the 



OO \ transition dynamics of RNA hairpin and increase the mean transition rate. Furthermore, the RNA 

dynamics can still discriminate a signal with resonance frequency even when the signal is mixed 
among other signals with nonresonant frequencies. In fact, our computational demonstration of 
thermally induced resonance in RNA hairpin dynamics is a direct realization of the phenomena 
called stochastic resonance (SR) and resonant activation (RA). Our study, amenable to experimen- 
tal tests using optical tweezers, is of great significance to the folding of biopolymers in vivo that 
are subject to the broad spectrum of cellular noises. 
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The cell is a highly crowded and dynamic place with millions of biomolecules and their 
interactions [l| . Only recently the notion that the self-assembly and function of biomolecules 
in vivo could be very different from those in isolation has been highlighted in terms of macro- 
molecular crowding and confinement [2 - 4 ] . However, the implication from dynamics in vivo 
should not necessarily be limited to the confinement or crowding effect. Besides the thermal 
noise, which is essential for activated internal motions and folding dynamics of biomolecules, 
of particular note are the nonequilibrium athermal signals produced from a plethora of ac- 
tive processes and collective actions among cellular constituents. For example, ribosome 
translocate along mRNA strand less than every 0.1 s per codon [ij]; Molecular motor, ki- 
nesin, takes step along the microtubule approximately every 10 ms jfj; Collective beating 
of cilia occurs approximately every 1 min 7]. How these nonequilibrium signals or noises 
affect fundamental processes in biology, particularly the folding dynamics of biopolymers 
8|, |9j, is the main concern of this paper. To address this issue, we study the effect of a weak 
oscillatory force on the transition dynamics of stretched RNA hairpins as a test case. 

In recent years a combination of experiments, theoretical arguments and simulations 
have been used to extensively study the conformational dynamics of biopolymers under 



mechanical control 



10M20| . Among others, paramount precision of optical tweezers, which 



can both exert piconewton force to molecules and probe structural changes in nm scale, not 
only allowed us to quantify thermodynamic and kinetic properties of the molecule at single 



molecule level 
of life 
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13| . but also provided unprecedented chances to manipulate the process 



24]. To emulate nonequilibrium athermal noise in vivo and study its effect on 
the conformational dynamics of biopolymers, high precision force instruments seem to hold 
good promise; yet, no experiment to date has explored this issue. 

Given that the dynamics of proteins and nucleic acids occur at overdamped fluctuating 
media, it may be surprising to learn that conformational transitions of biopolymers could 
become enhanced or synchronized to an externally exerted weak oscillatory signal. However, 
according to the notion of stochastic resonance (SR), sensitivity of a system to external 
signal can be improved under thermal fluctuation, so that the system can discriminate 



the weak information carrying signals 



25|. Theoretically, it has been shown for a single 



Brownian particle in a bistable potential that synchronization of barrier crossing events to 
an externally driven oscillation occurs not only at an optimal noise strength but also at 
resonance frequency 



25 



26(| . Furthermore, at another oscillation frequency that differs from 
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the frequency of synchronization, the mean first passage time over the barrier is minimized, 
which is called the resonant activation (RA) 2^]. Hence, it is possible to demonstrate the SR 
and RA in the conformational transition dynamics of biopolymers like RNA hairpins, whose 
folding landscape is bistable near transition mid-force, by properly tuning the frequency of 
additive oscillatory force. 



Since the pioneering work of Benzi et al. 28J, there have been efforts to extend the 



25 



301 ] . However, SR 



idea of SR to a vast range of physical and biological phenomena 
often mentioned in terms of macroscopic systems such as recurrence of the ice age cycle 
and sensory systems for higher order organisms, especially the latter of which is wired 
with complex biochemical networks, is always a subject of debate because the inherent 
complexity of a system consisting of a formidable number of variables often makes the origin 
of such resonance-like phenomena inconclusive. Although phenomenological modeling of 
system dynamics using SR is still possible, there seems to be little convincing proof that the 
biological systems actually use SR in the cellular level 31[. In this aspect, recent theoretical 
demonstrations of SR and RA of a polymer chain in a confined space has more direct 



relevance to the biology 32 



33f ] . Here we show that SR and RA can indeed be manifested 



in self-organization process of biopolymers by using molecular simulation. Although careful 
investigation is further required to confirm whether the SR and RA are actively utilized in 
the cell for its biological functions, it should be straightforward to demonstrate SR in the 
folding and unfolding dynamics of biopolymers by using single molecule force experiments. 
We provide a simple experimental protocol for the demonstration. 

RESULTS AND DISCUSSIONS 



Bistable transition of P5GA hairpin under constant tension : To study the 
dynamics of biopolymers under tension, we used the self-organized polymer (SOP) model, 
a new class of versatile coarse-grained model that have been proven powerful in describing 
diverse problems in biomolecular dynamics at single molecule level, which includes proteins, 



RNA, and molecular motors 



18-20 



34j-l36| . Among the molecules studied using SOP model, 



the 22-nucleotide P5GA RNA hairpin is the simplest one, whose dynamics was studied 
under various mechanical conditions such as force clamp and force ramp (See Methods), 
providing microscopic insights into RNA hairpins in single molecule force experiments 



3 



18H20|. Here, we performed Brownian dynamics simulation of the RNA hairpin under 
varying tensions (see Methods). As shown previously 3-20], the time trajectories that 



describe the dynamics of molecular extension (z) show that the RNA hairpin undergoes 
thermally-induced hopping transitions when the molecule is held around the transition 
mid-force f m . The force-dependent free energy profile F(z; f) = — ksT log P(z; f) with 
T = 300 K, calculated from the distribution of molecular extension P(z;f), becomes 
bistable at / = f m {F(z; f m ) = F m (z)) by forming the native basin of attraction (NBA) and 
unfolded basin of attraction (UBA) at zp ~ 2 nm and zjj ~ 8 nm, respectively (Fig. [Th.). 
The tension greater (or smaller) than f m , f > f m (or / < f m ) tilts the free energy profile 
towards the UBA (or NBA). The force-dependent probability of being in the UBA Pjj{f) 
is calculated using P v (f) = dze~ pF ^ / f£ dze~^ z ^ where /3 = (k B T)~\ z TS is the 
extension corresponding to the barrier top (transition state), z m and zm (see Fig. [TJl) are 
the minimum and maximum extensions, respectively. Pu(f), fitted to a two state model, 
P* = [i + e (/m-/)A z /fc s Tj-i gives f m _ 14 7 pN an ^ Az « 6 nm (Fig. [lb). From the time 

trajectories, the /-dependent mean folding (tf(/)) and unfolding time (ru(f)) are also cal- 
culated and found to "approximately" obey the Bell- like relation 37|] (Fig. [lb). It is of note 
that the deviation from linearity on a log folding time scale is due to moving transition state 



resulting from the underlying folding landscape that deforms under mechanical force 38|-|40|. 



Oscillation-induced coherent transition of the P5GA hairpin : According to 
theories, the nonlinear and bistable characteristics of the folding landscape of RNA hairpin 
guarantee the manifestation of SR and RA if an oscillatory driving is externally applied to 
the molecule just like the SR discussed in the context of Brownian particle on a bistable po- 
tential. However, given that there are many degrees of freedom whose dynamics experiences 
continuous dissipation over time, it is not obvious whether a small oscillatory force added 
to the RNA can indeed be sensed along the backbone of biopolymer and renders the folding 
transition synchronous to the oscillatory force. For proteins and RNAs where entropy play 
dominant role in shaping the underlying free energy landscape, the thermal noise strength 
(temperature), kept at a physiological condition, is not a good variable to vary. Instead, it 
would be more appropriate to vary the frequency of oscillatory force. To this end, we per- 
formed the many-body simulations using SOP model by applying a small oscillatory force 
5f sin (Qt) with 5f = 1.4 pN to the 5'-end of P5GA hairpin at / = 17 pN and monitor the 
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change in folding- unfolding dynamics, for various driving frequencies Q = 2tt/Tq. 

Our many-body simulations using the SOP model of RNA hairpin show that although 
the probability of being in the NBA at / = 17 pN is only « 5 % so that the folding 
is improbable, folding transitions occur more frequently and in coherence to a small 
oscillatory force at the period of Tq = T SR ~ 10.2 ms (Fig. |2b). Here, it is noteworthy 
that T SR for our RNA hairpin at / = 17 pN is similar to the sum of the two transi- 
tion times in the absence of oscillatory driving, 7>(/) + T u{f) ~ 10.8 ms (see Fig, [lb 
at / = 17 pN), a theoretically well-known SR condition for symmetric bistable potential [25]. 



Mapping P5GA dynamics onto one-dimensional free energy profile : To ex- 
plore the details of SR and RA in the dynamics of RNA hairpin more extensively, we 
simplified our many-body simulation using SOP model into one-dimensional (ID) simula- 
tion by showing that the dynamics of molecular extension z(t) from many-body dynamics 
can effectively be projected onto the dynamics of a quasi-particle under the ID free energy 
profile F(z;f). By using /-dependent free energy profile F(z; f) = F m (z) — (f — f m ) z 
where the F m (z) is directly interpolated from the energy profile at f m = 14.7 pN ob- 
tained by the SOP simulation (Fig. [TJ, we modeled an associated Langevin dynamics of 
the quasi-particle with a position z(t) as follows: jz = —dF(z; f)/dz + £(£). Here, 7 is 
the effective friction coefficient and £(t) is the random force that satisfies the fluctuation- 
dissipation theorem (£(£)£(£')) = 2kBT r y5(t — t') with (£(£)) = 0. Folding and unfolding 
times, determined from time trajectories z(t) on F(z;f), T F D (f) and t^- d (/) are in excel- 
lent agreements with those from the SOP model (tf(/) and Tu(f)) (Fig. [lb). In addition, 
by imposing reflecting boundary condition at z m ~ 1.01 nm and zm ~ 9.7 nm, which are 
the minimum and maximum molecular extension, we analytically calculated the mean first 
passage (Kramers) times Tp R (f) = dzeP F ^ z '^ f* M dz'e~^ F ^ z '''^ for the folding and 

T$ R (f) = J f£ dzeP F &fi j z zm dz'e-Wtf for the unfolding. We set the effective diffusion 
constant D = ksT/j w 13.7 /im 2 /s to match Tp R (f) (ru R (f)) (solid curves in Fig. db) to 
Tp(/) (ju(f)) at / = f m from the simulation results using SOP model (circles in Fig. [Tb). 

The excellent agreement among differently calculated mean transition times 7>(/) ~ 
T F D (f) ~ T F R (f) justifies mapping the many-body dynamics of P5GA hairpin under varying 
/ onto the dynamics on ID free energy profile described with molecular extension, z. In fact, 
previous studies showed that z is a good reaction coordinate on which both thermodynamics 
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and kinetics of P5GA hairpin are faithfully described as long as / does not greatly differ from 



fn 



19 



20] . This enables us to explore SR and RA of the RNA folding-unfolding dynamics 



in a simpler way without compromising the goal of our study. To take into account the 
effect of external oscillatory force, we simply added a term 5/sin(f2t); thus, the stochastic 
equation describing the dynamics of z(t) in the presence of external oscillatory signal reads 
as follows: 

BF (z) 

JZ = ^ + (/ - frn) + Sf Sin(m) + ^(t). (1) 



Resonant activation in the folding and unfolding transitions of RNA hairpins 

: To investigate the effect of oscillatory driving on transition kinetics, we calculated the 
average folding time Tp as a function of Tq for time trajectories generated using Eqfj] at 
/ = 17 pN for 5f = 1.4, 2.0 pN or unfolding time t v at / = 13 pN for 5f = 1.0, 1.4 pN 
(Fig. [3]). We found that the external oscillatory signal accelerates both folding and unfolding 
dynamics and maximizes the transition rates (defined by t f ,Ty ) at Tq = T RA rs 2.56 ms 
by about 60 % and 30 % for both folding and unfolding, compared with those in the absence 
of the oscillatory driving. The enhancement of folding (unfolding) rate could still be called 
RA, following the theoretical studies on symmetric bistable potential j^ . 

The presence of optimal condition of Tq = T RA that minimizes the transition time indi- 
cates that there should be a matching condition for two frequencies between conformational 
fluctuation in folding landscape, which is induced by thermal noise, and the oscillation 
of input signal. For Tq <C T^" 4 < r, net driving force provided by the fast oscillating 
force is essentially zero; ft +to Sf sin (2irt/Tn)dt ~ J °° T^Sf sin (27is)ds ~ 0. In contrast, 
when T RA <C Tq < r, the molecule experiences external tension that slowly varies be- 
tween f — Sf and / + 5f during the transition time r depending on the time interval 
t < t < t + r. As a result, for Tq w 100 ms, Tp at / = 17 pN and 5f = 1.4 pN 
varies between 2 ms and 30 ms (Fig. [lb), and the corresponding mean folding time becomes 
Tp = 7 Jt" +T Pfif (t))rp(f (t))dt ~ 9 ms as shown in Fig. [3^, which is greater than the 
t f w 8 ms at Tq = T RA [4l|. 

The extent of RA varies with the amplitude of oscillatory force (5f) (compare circles 
and squares in Fig. [3]). It is conceivable that 5f x 5z F where 5z F = z\j — zps ~ 4 nm at 



6 



/ = 17 pN (Fig. [T]) provides extra work to rectify the random fluctuation of molecule in 
UBA. Therefore, the ratio of 5f x 8z F to free energy barrier AF$ from the UBA to NBA, 
i.e., 5f x 5zp/AF$ should determine the enhancement of folding rate at the RA condition. 



Stochastic resonance in the refolding dynamics of RNA : Folding time distri- 
butions Pfoid(t) (Fig. [4b ) calculated from time trajectories (Fig. are useful measure for 
characterizing SR. At fast driving Tq = 0.16 ms (Fig. 0b-l), Pf id(t) decays exponentially 
with a decay time 10.5 ms, which is consistent with the folding time in the absence of 
driving, t f (see Fig. [Tb). But at Tq = 10.24 ms near the resonant condition, Pf id{t) dis- 
plays multiple peaks whose periodicity corresponds to the driving period of input signal, 
which is approximately equal to r F (Fig. Hb-2). The multiple peaks are superimposed on 
the above-mentioned exponentially decaying contribution (Fig. Hb-1) 25J. At a slow driving 
Tq = 40.96 ms (Fig. Hb-3), the Pf id{t) has a broad and small peak around 25 ms and shows 
only a weak signature of coherence. 

The degree of synchronization is quantified by calculating the P\ value corresponding to 
the area under the first peak of Pfoid(t), which is defined as P\ = JJ""™ 2 dtpf id(t) where the 
integration is done over a time interval between r m i ni and r m j n2 that defines the first peak in 



the Pfoidit) [25j, [26|, |42j. For example, for Tq = 0.16 ms, the Pfoidif) exponentially decreases 



with no peak, thus Pi = 0. In contrast, Pf id{t)s with Tq = 10.24 ms and 40.96 ms (Fig. Hb-2 
and Hb-3) have the first peaks around 7 ms and 25 ms, respectively, leading to non-zero P\. 
Pi(Pn) calculated for the data p/ id{t) of our ID simulation for P5GA hairpin (filled circles in 
Fig. Hb) is non-monotonic with Tq and is maximized around Tq = T SR « Tp + Ty « 10.8 ms. 
Notably, Pi(Tq) rise sharply for T n > T RA ~ 2.56 ms (Fig. 0b) and the T SR that maximizes 
Pi differs from the T RA that minimizes t>. The Pi for the unfolding with / = 13 pN and 
5f = 1.4 pN, by calculating unfolding time distribution p un foid(t), are also depicted with 
empty circles in Fig. HJi. We find that the unfolding is also synchronized with an optimal 
driving period around Tq ~ t f + ~ 8 ms (see Fig. [lb). Thus, RNA folding (unfolding) 
becomes maximally coherent at T SR ps t f + t^. Due to asymmetric shape of potential 
leading to Tp ^> tjj, the SR condition for folding becomes w r F + Ty ~ Tp, which is 
worth comparing with the conventional SR condition Tq t f + Ty = 2r F for symmetric 
bistable systems. 

There are many alternative measures to demonstrate the manifestation of 
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SR. One measure, directly amenable to experiment, is the power spectrum 
S(u) = f^ oo e~ t2nUT (z(t + r)z(t))dr where (...) denotes ensemble and long time aver- 
ages for the trajectory z(t) 25j. The S(u) shows a narrow sharp peak with the frequency 
v = 1/Tq ~ 1/10 ms" 1 (see the inset of Fig. Hb-2). In addition, spectral behavior of 
dynamic susceptibility defined x( u ) = ( z ( u )) I 'J "{ v ) > which also becomes optimal at v « 1/10 
ms -1 , can be used to quantify the response of system to an external driving force. 



Biopolymers can filter the resonance frequency from the spectrum of an 
external signal : Due to the multitude of cellular processes nonequilibrium noise in the 
cell is decomposed into various frequencies. To check if our model can discriminate the 
resonant frequency using SR, we performed the ID simulation for RNA hopping transitions 
under / = 17 pN by adding oscillatory forces consisting of different angular frequencies 
XliLi $fi sin(f2j£) where Q\ = 27r/0.1 ms" 1 , f2 2 = 27r/10 ms" 1 and f2 3 = 27r/100 ms" 1 ; here, 
the optimal folding period for the SR, Tsr ~ 10 ms, is included in the i — 2. The Pf id(t) for 
uniform amplitudes of the additive forces with Sfi = 6/2 = 6/3 — 1.4 pN (Fig. [5]) displays 
an oscillatory pattern with multiple peaks, suggesting that the folding is synchronous with 
the oscillatory force with the period T SR . Even for the case of nonuniform amplitudes, 
in which 5j\ and 6/3 have greater values than Sf2, the signature of synchronization to 
the oscillatory driving with T SR is retained (Fig. [5]). Furthermore, even when we choose 
5f randomly from the values distributed between and 2 pN at every time setp, we still 
found the Pf id(t) with multiple peaks, which suggests that the P5GA system is capable of 
filtering the signal with frequency of 2n/T SR out of the nonequilibrium noise. 



Experimental demonstration of SR in biomolecular dynamics : Although the 
results discussed above are based on the numerical simulation that applies mechanical force 
directly acting to the ends of RNA hairpin without mediating polymer handle as in actual 
optical tweezers experiment, our findings on SR and RA conditions for folding transition 
should not change qualitatively. As quantitatively discussed by Hyeon et al. 19[] as well as 



recognized by others 



111 ], the presence of handle polymers simply slows down the transition 



kinetics by pinning the motion of the two end points of biopolymers connected to the handles. 
Furthermore, in all practical situations, the conformational transition time of biopolymers 
is much slower than the tension propagation time 43(. Therefore, even in the presence of 
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handles SR is still realized at the optimal frequency, corresponding to the sum of the folding 
and unfolding transition times that are measured in the presence of handles. Many-body 
simulation using SOP model of RNA hairpin in the presence of 20 nm handles attached to 
the two ends of hairpin confirms the realization of SR in the presence of handles (see SI text 
and Fig.Sl). For the purpose of experimental demonstration of SR using biopolymers, it is 
recommended that the experiment is conducted at the transition mid- force (f m ), at which 



the folding and unfolding transitions are the most frequently observed [11|, [13] (Fig.S2). In 
this case, SR condition is expected to be formed around twice the mean hopping transition 
time (« 2r F (f = f m )). 



CONCLUDING REMARKS 



Originally, the SR refers to the phenomenological observation that a system under an 
optimal thermal noise can amplify a weak external signal if the system is nonlinear. This 
paper restates the SR and RA condition by focusing more on the transition dynamics of 
biopolymers rather than the amplification of signal. Furthermore, unlike the many stud- 
ies that consider the SR using noise strength, i.e., temperature, as a variable, the present 
study focuses on the "bona-fide SR" with respect to the driving frequency. For biologically 
relevant environment, whose temperature range should be narrowly defined, it is more ap- 
propriate to consider the effect of driving frequency on SR. Folding and unfolding transitions 
of biopolymers enhanced under and synchronized to an external periodic signals, which is 
clearly demonstrated in our simulations, suggest that as long as the SR condition is met, 
conformational sampling of biopolymers in vivo occurs more frequently than in vitro, in 
response to the spectrum of athermal noise in the cell. 

True manifestation of SR and RA in the context of biopolymer dynamics in vivo should 
hinge on existence of the actual periodic frequency that is in resonance to the transition 
dynamics of biopolymer. It would be of great interest to investigate the power spectra 
of cellular noise and how biological dynamics, for instance, conformational transitions of ri- 
boswitch, are modulated in the presence of athermal signals from other dynamic components 
such as motion of ribosome translocation and associated gene regulation fl^ |. 
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MATERIALS and METHODS 



Model : We used the SOP model to simulate the dynamics of P5GA RNA hairpin under 



tension. In SOP model, the effective Hamiltonian is defined as follows 

kRf 



H 



SOP 



5> 



N-3 N 

+ £!> 

i=l j=i+3 
N-3 N , 

i=l j=i+3 ^ 



(r i>i+1 - r° i i+1 ) 2 ' 



Rl 



hi 



12 



hi 



h3 



A 



N-2 
i=l 



r i,i+2 



(2) 



where we modeled each nucleotide as an interaction center; thus N = 22. The first term 
restrains the bond length between the neighboring beads by using the finite extensible 
nonlinear elastic potential with k = 20 kcal/(molxA 2 ), and Rq = 0.2 nm. Here, r^j+i is the 
distance between neighboring interaction centers % and i + and r° i+1 is the distance in the 
native structure. Lennard- Jones potential was employed to stabilize the base pair present 
in the native structure in the second term with = 0.7 kcal/mol, but make the non-native 
base pair repulsive in the third term, so that the model can easily reach the native topology 
without being trapped into non-native conformations. If % and j are in contact in the 
native state = 1; otherwise = 0. r°j is the distance between i and j interaction 
centers that satisfy \i — j\ > 2 and Ay = 1. For non-native pairs, the strength and range 
of interaction were set by e; = 1 kcal/mol and a = 0.7 nm. The last term involving the 
repulsion with a* = 0.35 nm maintains a certain angle between i and i+2 interaction centers. 



Simulation : We performed Brownian dynamics simulation by numerically integrating 
the discretized equation of motion 



r a>i (t + At) = f a>i {t) + 
n=i(t) = 



At ( dH 



SOP 



dr r 



(3) 



where a — x,y, z and the random force satisfies (£, a ,i{t)) = and it is assumed that 

the probability distribution of the noise in the random force is Gaussian, i.e., V[£ a ,i(t)] oc 



exp l — Jdr 



and (£ a:i (t) ■ £pj(t + n A*)) = l&fS^JiJa^ with 
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0,1,2,.... We 



10 



probed the folding- unfolding dynamics of RNA hairpin by applying tension / to the 3 '-end 
of the hairpin while fixing the position of 5'-end at temperature T = 300 K. Also, we chose 
the friction coefficient £ = lOOm/r^ with tj, = (ma 2 / e/J 1 / 2 ~ 3 — 5 ps corresponding to the 
value that a nucleotide whose size a ~ 0.5 nm feels in water environment whose viscosity is 
~ 1 cP [181 ] . We applied the tension / between 13 < / < 18 pN (Fig. to the molecule 
and added small sinusoidal forces of Sf sin (fit) with varying fl values (Fig. [2]). We set 
the simulation time step At = 4.7 x 10~ 5 /isec and recorded data at every 10 5 time steps 
(4.7 /isec). After equilibrating the RNA in native state for ~ 47 fis, we applied tension / or 
/ + Sf sin (fit) and recorded the z(t) of each trajectory for t > 10 2 ms. 
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FIG. 1: Thermodynamics and kinetics of P5GA hairpin under tension, (a) Free energy profile 
as a function of molecular extension z under tension /, obtained from the molecular simulation 
using SOP model. Shown are the three free energy profiles at / = 13.0, 14.7, and 17.0 pN. At 
/ = fm = 14-7 pN, the probability of being in UBA and NBA are equal. The profile is tilted towards 
UBA (or NBA) when / = 17.0 > f m (or / = 13.0 < f m ). (b) Unfolding and folding probability 
Pu(f) and Pf(/)(= l-Pu(f)) of P5GA against tension /. Results from the SOP model (P F (f) and 
Pjj{f)) and two state model (solid curves) are in good agreements. P^ D (/) (or Pp D {f)) calculated 
from the tilted ID free energy profile interpolated from results of the SOP model displays a good 
agreement with other results, (c) Mean folding time (tf(/)) and unfolding time (rf/(/)) at varying 
tensions, obtained from the SOP simulation (circles), T F D {f) and rlj D {f) obtained from the ID 
simulation on the tilted free energy profile (squares), and mean first passage times (rp R , T^j R ) 
calculated using the tilted free energy profile (solid curves) are in a good agreement. The transition 
mid- force determined at the intersection between rp(f) and Tu(f) is f m = 14.7 pN. 
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FIG. 2: Demonstration of SR using simulation of P5GA hairpin under an oscillating force, (a) 
Structural ensemble of P5GA hairpin at NBA and UBA shown on the free energy profile F(z; f = 17 
pN). Time trajectories of molecular extension z(t) at constant tension / = 17 pN are shown 
from the simulation of SOP model when a small time dependent oscillatory force 5f sin (27rt/Th) 
(5f = 1.4 pN) is added with (b) r n (= 2vr/0) = 1.3 ms, (c) 10.2 ms, and (d) 82 ms. Coherent 
hopping transition between folding and unfolding is found with Tn = T SR ~ 10.2 ms, which 
clearly shows the SR. The oscillatory forces with Tq = 10.2 ms and 82 ms are overlaid on the time 
trajectory in c and d, respectively, to emphasize the effect of synchronization. See Fig. SI for the 
SR of P5GA hairpin in the presence of handles. 
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FIG. 3: Resonant activation (RA) in folding and unfolding transitions of P5GA. (a) Mean folding 
time Tp of P5GA under / = 17 pN as a function of Tq with 5f = 1.4 pN (filled squares) and 
5f = 2.0 pN (filled circles). The transition times are minimized at Tq = T RA ~ 2.56 ms. (b) Mean 
unfolding time t\j of P5GA under / = 13 pN as a function of Tq with 5f = 1.0 pN (empty squares) 
and Sf = 1.4 pN (empty circles). The transition times are minimized at Tq = T RA ~ 2.56 ms. 
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FIG. 4: SR in the folding and unfolding transitions of P5GA under time dependent tension 
/ + 6fsin(p,t). (a) Time trajectories and (b) folding time distributions Pf id(t) with / = 17 pN 
and 5f = 1.4 pN at Tq = 0.16 ms, 10.24 ms, and 40.96 ms are shown. The synchronization (SR 
condition) occurs around T n = T SR « t° f + r£ » 10.8 ms, as shown in (a-2). The inset in (b-2) 
depicts the power spectrum SfV), showing the sharpest peak at v = 1/T SR ~ 0.1 ms -1 . (c) Mean 
folding time rp and measure of coherence Pi as a function of Tq. (d) Mean unfolding time tjj and 
measure of coherence -Pi as a function of with / = 13 pN and Sf = 1.4 pN. Note that the SR 
and RA conditions are made at different Tqs in (c) and (d). 
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FIG. 5: Filtering of SR condition. Shown are the Pfoid(t) of P5GA under / + Y2i=i $fi sin(f2jt) at 
varying combinations of 8fi with = 27r/0.1 ms~ , — 27r/10 ms -1 , = 27r/100 ms -1 and 
/ = 17 pN. For filled symbols, we fix 8/2 = 1.4 pN, corresponding to the amplitude for the optimal 
frequency ^2 for SR, and vary the other amplitudes to be 8fi = Sfy = 1.4 pN (filled squares), 
1.7 pN (filled circles) and 2.0 pN (filled triangles), respectively. The starred symbols represent the 
cases where 8f±, 8f% and 8fs are taken randomly between and 2 pN every time step to make 
a nonequilibrium noise 8f{t). Notably, the folding transition filters the optimal driving period 
despite other drivings with larger or transient amplitudes. 
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SUPPORTING INFORMATION (SI) 

Stochastic resonance in the presence of handles : To study whether the SR is still 

realized when the weak signal is indirectly exerted to the system through handle polymers, 
just like in the real experimental condition in laser optical tweezers (LOTs), we simulate 
folding dynamics of an RNA hairpin by attaching stiff handles to the two of the P5GA RNA 
hairpin modeled with SOP representation (Eq.[2] in the main text). Each handle is modeled 
by using the following energy potential [lj]: 



The first term denotes the chain connectivity, with an equilibrium distance r$ = 0.5 nm 
and spring constant k$ = 1.4 x 10 4 pN ■ nm -1 . The second term is the bending penalty 
potential with bending modulus k^ = 561 pN-nm, which gives rise to the persistence length 
l p = 70 nm. We take N = 40 for the length of handle, which amounts to the contour length 
of L = 20 nm. Fig. [Si] shows time trajectories of the extension of the entire system, which 
includes the extension of two handles, at / = 17 pN and 5f = 1.4 pN, for three different 
driving periods. The molecular extension z(t) (left) and the total extension zi(t) (right) 
are depicted in the red and black axes, respectively. As shown in the trajectories, the phase 
of Zi(t) is compatible with that of z(t). Even in the presence of handles, the SR is realized 
at an optimal period 60 ms (Fig. ISTb) where the folding transitions are synchronized to the 
frequency of external driving. The other trajectories under relatively fast (Fig. IS lk ) and 
slow (Fig. ISTb ) driving show random transitions regardless of the drivings. As it is discussed 
in the main text and in the Ref. [lj], the presence of the handles pins the overall molecular 
motion and slows down the transition kinetics of the molecule sandwiched between the 
handles. We have estimated that the transition dynamics with handles with L = 70 nm and 

n 

L = 20 nm tend to be almost about six-fold slower pj than the Kramers times (~ 10 ms for 
folding and ~ 0.7 ms for unfolding) without handles which we have considered in the text, 
and found SR (Fig. ISTb) at an optimal driving period Tq rs 60 ms which is approximately 
sum of the Kramers times. 

Stochatstic resonance at the the transition mid-force : To demonstrate the SR using 
LOTs experiment, it would be most straightforward to perform the experiment around the 
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external force value near the transition mid-force, at which the transition from one basin 
to another is most frequent. To show SR under the transition mid-force (f m = 14.7 pN), 
we have performed the SOP simulations for hopping (folding and unfolding) dynamics of 
a P5GA under f m . Fig. [S2] shows time trajectories of the P5GA extension z(t) for three 
different driving periods, at 5f = 1.4 pN. SR is clearly shown in the Fig. IS2b at an optimal 
driving period Tq = 4 ms. In this case, the optimal driving period is found to be around 
twice the mean hopping transition timee (~ 2t>(/ = f m ) ~ 4 ms, see also Fig.lc where 



t~f = t~u =2 ms), in consistent with t 
been discussed in the SR community 



re time-scale matching condition as has extensively 
2j-|4j. The other trajectories, under relatively faster 
(Fig. IS2~k ) and slower (Fig. IB2~b ) drivings which deviate from the optimal frequency, show 
random transitions. 
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FIG. SI: Demonstration of RNA hairpin trajectories under SR condition in the presence of handles. 
Together with the SOP representation of RNA hairpin, we considered the handle of contour length 
L = 20 nm and persistence length l p = 70 nm. The time trajectories of molecular extension z(t) 
(red;left axis) and the extension with handles zi,(t) (black;right axis) are displayed at constant 
tension / = 17 pN, driving force 5fsm(£lt) with 5f = 1.4 pN, and period Th = 2n/Cl. The phase 
of extensions Zl(£) overlaid on the z(t) is in good agreements. Both z(t) and zi,(t) show random 
transitions at (a) = 10.0 ms and (c) Tq = 150.0 ms. (b) The folding transitions are synchronized 
with an optimal driving period Tq = 60.0 ms. The blue curves are overlaid to show the phase of 
oscillatory driving forces. 
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FIG. S2: Demonstration of bona-fide SR using the SOP model under the transition mid-force 
f m = 14.7 pN in the absence of handles. We depicted time trajectories of molecular extension z(t) 
at f m and <5/sin(f2i) with Sf = 1.4 pN and period Tq = 2-rr/O,. z(t) show random transitions at 
(a) — 0.1 ms and (c) = 10.0 ms, while coherent transitions occur at an optimal period (b) 
Tq = 4.0 ms. The blue curves are overlaid to show the phase of oscillatory driving forces. 
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